Dynamics of Scalar Field Dark Matter With a Cosh-like Potential 



On 
O 
O 
(N 

(N 

o 

u 

Or 
i < 

o ■ 



> 

o\ 
m 
o 

\6 
o 

OS 

o 



X 



Tonatiuh Matos, 1 '^ Jose-Ruben Luevano, 2, Israel Quiros,^R 
L. Arturo Urena-L6pez*, 4, and Jose Alberto Vazquez* 1 '^ 

1 Departamento de Fisica, Centra de Investigation y de Estudios Avanzados del IPN, A. P. 14-740, 07000 Mexico D.F. 

Sj \ Departamento de Ciencias Bdsicas, UAM-A, CP. 02200 Mexico, D.F., Mexico. 
Universidad 



Mexico. 



Central de Las Villas, Santa Clara, CP 54830, Cuba 
Departamento de Fisica, DCI, Campus Leon, Universidad de Guanajuato, CP 37150, Leon, Guanajuato, Mexico. 

(Dated: June 2, 2009) 

The dynamics of a cosmological model fueled by scalar field dark matter with a cosh-like potential 
plus a cosmological constant is investigated in detail. It is revealed that the late-time attractor is 
always the de Sitter solution, and that, depending on the values of the free parameters, the oscillating 
solution of the scalar field - modeling cold dark matter - mediates between some early stage (say, 
the radiation-dominated solution) and the accelerating de Sitter attractor. 
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I. INTRODUCTION 

One of the greatest mysteries of modern Cosmology 
doubtless is the nature of dark matter (DM) (for a review 
seep]]). We know that more than 23% of the matter of the 
universe is attractive, but of unknown nature. The most 
accepted candidate are particles from the Minimal Super- 
symmetric Standard Model (MSSM)0. This paradigm 
has been proved at cosmological level with a great suc- 
cess. For example, it predicts very well the formation and 
clustering of galaxies and the micro-Kelvin temper- 
ature fluctuations of the cosmic microwave background 
(CMB) of the universe[ij]. These two predictions repre- 
sent the main achievements of the model. 

Nevertheless, in the last decade, some inconsistencies 
of the model when calculations are compared with the 
observations at galactic level, have arosen. For example, 
amongst other problems, the predicted number of satel- 
lite galaxies around big galaxies is much bigger than the 
observed oneQ, and the DM density profile at the center 
of galaxies seems to be steeper than observed @. 

Even though there have been several attempts to deal 
with these inconsistencies Q, one compelling possibility 
is the hypothesis that the DM is a scalar field, more than 
a Standard Model (SM) particle The idea is rooted 
on the fact that almost all the unified theories of physics 
contain scalar fields as the simplest geometrical objects. 
These scalar fields are often called Higgs, Infiatons, Dila- 
tons, Scalerons, Radions, etc. All of them are needed to 
achieve mathematical and physical consistency of these 
theories. With the discovery of the dark energy (DE), at- 
tempts to relate scalar fields with the DE also have had 
some success (see@ for a comprehensive review). 
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In Ref. [8| some of us proposed that a scalar field rolling 
down a convex self-interaction potential can be a reason- 
able candidate for DM, calling this hypothesis as Scalar 
Field Dark Matter (SFDM). This hypothesis has some 
nice features. For example, the SFDM does not need ex- 
tra hypothesis to explain the flat DM profile in the center 
of galaxies [l(| , or the number of satellite galaxies around 
the Milky Way [HI- 

The SFDM hypothesis has been investigated for a 
number of self-interaction potentials like the cosh-likc 

nana 

ITU ], and the quadratic ones [3, [l5j], with the 
consequent discovery of several exact solutions of cosmo- 
logical interest. However, within the cosmological con- 
text, a full and detailed study of the dynamics of SFDM 
models, to uncover their relevant asymptotic properties, 
is still desirable. 

The main goal of the present paper is, precisely, to 
study within the cosmological context and by means of 
the dynamical systems tools, the asymptotic properties of 
the SFDM model driven by a cosh-like potential. Several 
relevant cosmological solutions will be correlated with 
such important dynamical systems concepts like past and 
future attractors, signaling the way the cosmic dynamics 
transits from early-time to intermediate and then to late- 
time asymptotic states. 

It will be demonstrated, in particular, that the SFDM 
model driven by a cosh-like potential with a cosmologi- 
cal constant term added, is a good candidate to correctly 
describe a cosmic dynamics whose fate is, eventually, the 
inflationary de Sitter regime that can be correlated with 
present accelerated stage of the expansion of the Uni- 
verse. Due attention will be paid to the late-time oscilla- 
tory solution that can e associated with cold dark mat- 
ter (CDM) behaviour in the model, and that, inevitably, 
leads to the inflationary de Sitter attractor. 

The paper has been organized as it follows. The rele- 
vant physical features of the model are discussed in the 
next section, and then, in section III, its mathematical 
features are specified. The details of the (linear) dynam- 
ical systems study of the SFDM model with a cosh-like 
potential, are given in section IV. The main results of 
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the study are summarized in section V. Finally, in sec- 
tion VII, the physical discussion of the results and brief 
conclusions of the study are provided. We use the natural 
system of units where k 2 = 8irG = c = 1 . 



II. SCALAR FIELD MODEL WITH A 
COSH-LIKE POTENTIAL 

In a cosmological context, scalar fields with a cosh- like 
self-interaction potential have been studied, for instance, 
in[13]. In that reference the following potential was in- 
vestigated: 



V(4>) = F o [cosh(A0) - 1] / 



(1) 



where A and (3 are free parameters. It was shown, see 
alsofHj], that during the oscillatory phase the virial theo- 
rem gives the following expression for the mean equation 
of state: 




{3-1 
13+1 



(2) 



In consequence, for (3 — 1 the scalar field behaves like 
pressureless dust = 0. A scalar field potential with 
this value of (3 could therefore play the role of cold dark 
matter in the universe. For j3 < 1/2 this potential is a 
good candidate for the quintessence instead p^j. 

In this paper we will focus our attention, precisely, in 
a cosh-like potentialfl]) with (3 = 1 with a cosmological 
constant term added to it. The latter ingredient is nec- 
cessary to include dark energy in the model, as long as os- 
cillations of the scalar field around the minimum of cosh- 
potential play the role of the CDM (in fact SFCDM). 

A reasonable physical assumption for the initial con- 
ditions is that the energy density in the scalar field was 
comparable to that of radiation at very early times, say, 
at the end of inflation [l3[ ■ At such early times, the large 
value of <p makes the cosh-like potential^ behave as an 
exponential one of the form V(cj>) oc exp(— f3\(j>), and the 
scalar field energy density decreases rapidly so that it 
remains subdominant until recently. 

In later times the form of V(<p) changes to a power law, 
resulting in rapid oscillations of the scalar ield around 
(j) = 0. At this stage the scalar field equation of state 
mimics that of CDM in agreement with the present cos- 
mological standard model paradigm. Due to the addition 
of the cosmological constant term, the resulting picture 
does indeed represent a unification of the dark matter 
and the dark energy. The present model might be called 
A-SFCDM. 

The goal of this paper will be to describe the dynamical 
picture produced by the cosh-like potential {J) with (3 = 
1 and with the addition of a cosmological constant term, 
to put the above physical discussion on solid dynamical 
systems grounds. 



III. MATHEMATICAL FEATURES OF THE 
MODEL 

In this paper we focus in a Friedmann-Robertson- 
Walker space-time with flat spatial sections, filled with 
a mixture of two fluids: i) a perfect fluid of ordinary 
matter with density p 7 and barotropic index < 7 < 2 
(7 = 1 for dust, 7 = 4/3 for radiation, etc.), and ii) a 
united A-SFCDM component with energy density — 
(p 2 /2 + V(<p) and parametric pressure p^, = (p 2 /2 — V((p). 
The relevant cosmological equations of the model are the 
following 



H = 



(p +7y9 7 



-3*0-^, 



p 7 = -3i?7p 7 , 
plus the Friedmann constraint: 

If 1 :< 



H 



Pi 



V 



(3a) 

(3b) 
(3c) 

(4) 



As anticipated, here we deal with a cosh-like potential 
of the form 



V(cp) = Vb[cosh(A0) - 1] + A, 



(5) 



where A > is a positive cosmological constant, and the 
parameter A is also a possitive quantity (A > 0). As it has 
been discussed in the former section, the potential V((j>) 
in ([5]) comprises both CDM and dark energy in a united 
ansatz. Actually, for a vanishing A = 0, this potential is 
a particular case of JT]), when the free parameter [3 = 1, 
meaning that, once the oscillatory phase is attained, the 
scalar field driven by the potential ([5|) behaves like dust, 
i. e., just like CDM (see equation @). Since in this 
model A ^ 0, the scalar field energy density performs 
damped oscillations around V m i n — A, meaning that the 
CDM energy density - accounted for by the oscillatory 
component - decreases until, eventually, the cosmological 
constant component dominates. 

The potential V{4>) in © is a minimum at <fi = 0, =>■ 
Vmin = V(Q) = A. Additionally, in a natural scenario for 
the cosmic dynamics driven by V((p), the scalar field cp 
runs from arbitrarily large negative values (\<p\ ^> 1), to 
vanishing ones (\(p\ 1). In consequence, at early times 
the dynamics is driven by an exponential potential 



H » l/A 



V{cp) a V e 



-\<t> 



V 



V 



(6) 



whereas at late times it is associated with a power-law 
potential plus a cosmological constant: 



1^1 « l/A 



V(4>) a ^ 



A. 



V X 2 . (7) 



The latter is just the quadratic potential that has been 
formerly studied, for instance, in references [3, [HI]. 



3 



IV. DYNAMICAL SYSTEMS STUDY 

The dynamical systems tools offer a very useful ap- 
proach to the study of the asymptotic properties of the 
cosmological models 0, [l7| ■ In order to be able to ap- 
ply these tools one has to (unavoidably) follow the steps 
enumerated below. 

• First: to identify the phase space variables that 
allow writing the system of cosmological equations 
in the form of an autonomus system of ordinary 
differential equations (ODE). There can be several 
different possible choices, however, not all of them 
allow for the minimum possible dimensionality of 
the phase space. 

• Next: with the help of the chosen phase space vari- 
ables, to build an autonomous system of ODE out 
of the original system of cosmological equations. 

• Finally (some times a forgotten or under- 
appreciated step): to indentify the phase space 
spanned by the chosen variables, that is relevant 
to the cosmological model under study. 

After this recipe one is ready to apply the standard tools 
of the (linear) dynamical systems analysis. 



A. Autonomous System of ODE 

Let us to introduce the following dimensionless phase 
space variables in order to build an autonomous sys- 
tem out of the system of cosmological equations (|3c|) 
and ®[18] 



, y 



V 



V3H 



(8) 



After this choice of phase space variables we can write 
the following autonomous system of ordinary differential 
equations 

X ' = -Vf^y 2 -3z+~z(2x 2 + 7 Q 7 ), (9a) 



, /3 <9aV 3 ,„ , „ , 

y = \J 2^r x v + 2 v[2x + ^ ' 



(9b) 



where a prime denotes derivative with respect to the 
time variable r = In a (properly speaking the number of 
c- foldings of expansion), and the dimensionless density 
parameter f2 7 = p~ ( /3H 2 is given through the following 
expression 



n 7 = 1 - x 2 - y 2 



(10) 



which is just a rewriting of the Fricdmann constraint l|4|). 

As long as one considers just constant and exponential 
self-interaction potentials (d^V — and d^V — const, re- 
spectively), Eqs. (|9"aj) and (|9L>[) form a closed autonomous 



system of ODE. However, if one wants to go further and 
to consider a wider class of self-interaction potentials be- 
yond the exponential one - as it is the case in the present 
study -, the system of ODE (|9a[) and (|9b[) is not a closed 
system of equations any more, since, in general, d^V is a 
function of the scalar field <f>. 

A way out of this difficulty can be based on the method 
developed in (l9| . In order to be able to study arbitrary 
self-interaction potentials one needs to consider one more 
variable v, that is related with the derivative of the self- 
interaction potential through the following expression 



;A = -d^V/V 



,lnV. 



(11) 



Hence, an extra equation 

v' = -V6Xxv 2 {T- 1), (12) 

has to be added to the above autonomous system of equa- 
tions. The quantity T = V djV / (d<f,V) 2 in Eq. (p)) is. in 
general, a function of (f>. The idea behind the method 
m [H is that T can be written as a function of the vari- 
able v, and, perhaps, of several constant parameters. In- 
deed, for a wide class of potentials the above requirement 

- T = I» -, is fulfilled. 

Let us introduce a new function f(v) — v 2 (T(v) — 1) so 
that Eq. (|12[) can be written in the more compact form 

v' = -V6Xxf(v) . (13) 

Eqs. (|9ap. (|9bp. (jTDJ) , and (fT3"|) form a three-dimensional 

- closed - autonomous system of ODE 



Xy 2 v-3x + -x(2x 2 (14a) 



3 3 

-Xxyv + -y(2x 2 +^n j ) , 



= — V6Xxf(v) , 57 



1 



y 



(14b) 
(14c) 



that can be safely studied with the help of the standard 
dynamical systems tools [lTj. Notice the obvious symme- 
try of the ODE's (|14c|) under the simultaneous change of 
sign of A and v. X — > —A, v — * —v. 

An important limitation of the approach explained 
above is related with the fact that, for potentials that 
vanish at the minimum - such as, for instance V — 
V" o (cosh(A0) - 1)0, or V = m 2 <j) 2 /2 -, the variable 
v = ~d < pV/(XV) is undefined at this important point, 
usually important for the late-time dynamics. However, 
thanks to the non-vanishing cosmological constant term 
A in Eq.©, the latter is not the case in the present study, 
so that the approach of[19( can be safely applied. 

In terms of the phase space variables x, y, and v, the 
following relevant magnitudes: i) the deceleration param- 
eter q = — (1 + H/H 2 ), and ii) the effective equation of 
state parameter of the scalar field 74, = 2</> 2 /(</> 2 + 2V), 
can be written, respectively, as 



q = -1 



> 2 



2x 2 



(15a) 
(15b) 
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FIG. 1: A plot of function v(<j>) for the chosen values of the 
parameters of the potential A = 1, a = 10 - darker curve, 
a = 1 - dark-to-gray curve, and a = 0.5 - soft-gray curve, 
respectively. 
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FIG. 2: A plot of function cosh(A^)) vs v, for the chosen values 
of the free parameters of the potential (A = l,a = 10) - left- 
hand panel, and (A = l,a = 0.5) - right-hand panel. Both 
branches of the cosh-function, the positive (dark curve) and 
the negative one (gray line) are shown. Notice that, for the 
second case (a < 1), the whole cosh-function behaviour is 
depicted by a union of part of the negative and the entire 
positive branches. 



B. Function f(v) for the cosh-like Potential (T5j) 

We now derive the function f(v) in Eq. (fTS"]) (also 
present in Eq. (|14c|) ). for the potential of interest here 
(Eq.©). According to the definition (fTTj) for the cosh- 
like potential, one has 



v(cf>) 



sinh(A0) 



(16) 



cosh(A0) — 1 + a ' 
where a = A/Vq. From (|16p it follows, in particular, that 
lim v{(j>) = ±1, lim v((f>) = . (17) 

(p — >±oo tf> — >0 

The value a = 1 is a critical one. Actually, the func- 
tion (1161) is an extremum whenever: 



cosh(A0) 



1 



1 — a ' 



a(2 - a) 



(18) 



The fact that the cosh function is always equal or big- 
ger than unity means that v(4>) is an extremum only if 
< a < 1. For a > 1 the function v is, at least, a mono- 
tonically non-decreasing function as the scalar field takes 
values in the interval — oo < cf> < oo, therefore, in this 
case v e [—1,1]. 

In Fig. Q] the function v(4>) is plotted for the following 
sets of values of the parameters of the potential V(<p): 
(A = l,a = 10) - dark curve, (A = 1, or = 1) - dark- 
to-gray curve, and (A = l,a = 0.5) - soft-gray curve, 
respectively. As seen, in the latter case, since a < 1, 
v((f>) is a minimum at 

X(f) = — arccosh( ) = — arccosh(2) , (19) 

1 — a 



where v m i n = — 1 / \J ot{2 — a) = — -\/4/3, while at X(f> 
arc cosh(2) it is a maximum instead: v max — \J 4/3. 
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FIG. 3: A plot of function / + in Eq. (|22|l vs v, for the chosen 
values of the free parameters of the potential in Figfl]- left- 
hand panel. In the right-hand panel we show a plot of / vs 
v for a = 0.5 < 1 (A = 1). Both branches of f{v) in Eq.j22j, 
the positive (darker curve) and the negative one (gray line), 
are depicted. 



The function v(4>) in Eq. (|16p can be inverted to have 
the cosh function written in terms of 



(a - l)v 2 ± a/T+ a(a - 2)v 2 , . 

coshA^u = " " r 2 — > 20 

1 — v z 

where the "±" signs refer to the two different branches 
of cosh(A0)(t>). 

For values a > 1, the mathematical requirement that 
the cosh-function has to be, neccessarily, bigger than or 
equal to the unity, leads to the following condition 



± v/l + a(a-2)w 2 > 1, => |v| < 1 , 



(21) 



so that, as already mentioned before, only values of v 
obeying \v\ < 1 can be considered. This point is rel- 
evat when discussing the definition of the phase space 
since v is one of the variables spanning the phase space 
of the model under consideration here. For < a < 1 
the reasoning is not so simple, however, it is clear from 
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the former analysis that, in this case, \v\ < \v ext \ = 
l/y/a(2 - a). 

In order to illustrate the discussion above, in Fig. [2] 
we plotted the cosh-function vs v for two sets of values 
of the parameters of the potential: i) (A = 1, a = 10) 
- left-hand panel, and ii) (A = 1, a = 0.5) - right-hand 
panel, respectively. In both cases the darker curve is 
for the positive branch ("+" sign in Eq. (|20p ). while the 
gray line is for the negative branch ("-" sign in Eq. ([2H]) ) 
instead. 

Notice from the left-hand panel of Fig. [2] (a = 10 > 1), 
that only the positive branch actually depicts the right 
(whole) cosh-function behaviour. In the right-hand panel 
(a = 0.5 < 1), instead, the whole cosh-function be- 
haviour can be covered by a union of that part of the neg- 
ative branch to the left of the vertical asymptote v = — 1 , 
starting at infinitely large values of cosh(A0) (i. e., in- 
finitely large negative values of <p) until the gray curve 
meets the dark one (the positive branch) at the union 
point at v = v min = -l/ x /a(2 - a). 

The same is true for positive values of </>: the right 
cosh-function behaviour is depicted by a union of that 
part of the negative branch to the right of the vertical 
asymptote at v = 1, starting at infinitely large posi- 
tive values of <p, until it joints the positive branch at 
v = v max = 1/ y/a{2 — a). The range of intermediate- 
to-small values of <f) (including the beginning at <j> = 0, 
where the potential V(<j>) in Eq. ([5]) is a minimum), is 
completely covered by the positive branch of the cosh- 
function ([20]) . 

For the potential ((5|) the function f(v) appearing in 
Eqs. (fTB"]) and (|14c[) . can be written in the following way 



/» = 



± 1 



(1 - v 2 ) y/1 + a(a - 2)v 2 



(22) 



., 1 ± a/1 + a(a - 2)v 2 

where the ± signs refer to the positive and the negative 
branches of the function f(y) (directly related with the 
positive and the negative branches of the cosh-function 
in Eq. P0)) ). respectively. In Eq. (|2"2"|) the choice of the 
"+" or the "-" sign is to be made simultaneously in the 
numerator and in the denominator on the right-hand- 
side. Hence, for instance, for the positive branch of f(v) 
one has 

/>} _ + "--; a »^ , ete . ,23) 



a — 1 



It has to be emphasized that, while for values of the free 
parameter a > 1 the positive branch f + (v) in Eq. ([22]) 
is enough to depict the whole dynamics, for < a < 1, 
instead, one needs to contemplate both branches f^iv). 

Actually, in the later case the piece of the dynamics in 
the w-interval 

1 < H < , = V ( 24 ) 
y?a(2 - a) 

is uncovered by the correponding values of f~{v) in 
Eq. (f2"2")l : > 1). The rest of the dynamics - in- 

cluding the late-time behaviour - is uncovered by the 



piece of the positive branch liying in the interval \v\ < 1: 
/ + (M < !)• In consequence, for < a < 1, the whole 
dynamics is uncovered by 



fo 



<a<l 



r(M>i)u/+(M<i). 



(25) 



To summarize, for a > 1 the cosmic dynamics driven 
by the potential ([5]) can be associated with the follow- 
ing 3-dimensional (compact) phase space, spanned by the 
variables x, y, and v (we take into account only expand- 
ing cosmologies H > 0) 

* a >i = {(x,y,v) : < x 2 + y 2 < 1, 

M <1, 0<y< 1, |v| <1}, (26) 

and we have to worry only about the positive branch 
of f(v) in Eq.([22]). Meanwhile, for < a < 1, the 3- 
dimensional (compact) phase space is given by 



<Q<1 



{(x,y,v) :0<x 2 + y 2 < 1, 
|a?| < 1, 0<y < 1, |«| < 



y/a(2-a) 



W) 



and, depending on the piece of the dynamics one is inter- 
ested in, one has to rely either on the negative branch of 
Eq. (|2"2"l) /~(|w| > 1) (early times-to-intermediate dynam- 
ics), or on the positive brach instead / + (M < 1) (late- 
time dynamics), the whole dynamics being described by 
/o< Q <i in Eq.pSJ). 

In what follows, depending on the value of parame- 
ter a = A/Vb, we shall look for fixed points of the au- 
tonomous system of ODE (|14c[) . with f(v) given, either 
by: i) the positive branch of equation (|22[) . within the 
phase space ^ a >i defined in Eq. (f2"6")l . whenever a > 1, 
or ii) by fo< a <i defined in Eq. (|25[) . within the phase 
space ^a <a <i defined in Eq. (|2"T|) , if < a < 1. 



C. Equilibrium Points and Stability 

The fixed points of the autonomous system of 
ODE (|14c|) . in the phase space <3> a >i defined by Eq. ((26)) 
or in ^o<a<i defined by Eq. (|2T|) . are listed in Table Q] 
while the eigenvalues of the corresponding linearization 
matrices are shown in Table [TTJ 

In the case when the parameter < a < 1, by look- 
ing at the definition of the function f(v) in Eq. (|22[) . 
it might seem that, in addition to the critical points in 
Table HI there can be also equilibrium points associated 
with the values v = ±1/ \Ja(2 — a). However, if one 
looks at Fig. [3] (right-hand panel), one can see that at 



v = ±1/ y/a(2 — a) the derivative of the function /o< Q <i 
in Eq. (|25|) . is undefined, so that the linear approach un- 
dertaken in this investigation is not valid any more. This 
means that the above points in phase space might not be 
actual critical points, so that we shall not include them 
in our analysis. 

Existence of the matter-dominated solution (equilib- 
rium point Pi in Table |TJ , is independent on the value of 
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TABLE I: Properties of the fixed points for the autonomous system (|14c 
Existence fi 7 fl^ 7^ 

P x v Always 1 undefined -1 + ^ 

P 2 1 1" 12 2 



P3 
P 4 



P 5 10 Always 10 -1 



f2 


O 


1 

1 


U 





1 





1 




33 




A* 





1 



TABLE II: Eigenvalues of the linearization matrices corresponding to the fixed points in table [I] We have used the following 
definition R = ^/lhy 2 - 287 + 12 (1 + 2 7 2 [37 - 2 + 12(1 - 7)7 2 /A 2 -7 2 (1 - 3-y 2 /X 2 )] /A 2 ). 



Pi Ai A2 A3 



Pi 


-3 + 37/2 




3 7 /2 





Pi 


6-3 7 




3- \/372A 


VE\ 


p 3 


-3 + A 2 /2 




A 2 -37 


A 2 


Pi 


-3(2 - 7 )/4 + V3-R/4 


-3(2 


- 7 )/4 - VSR/4 


3 7 


Pr, 


-37 


-3(1 


- v/1 - 4A 2 /a)/2 


-3(1 + y/1 - 4A 2 /a)/2 



the variable v, meaning that this phase of the cosmic evo- 
lution may arise at early-to-intermediate times, as well 
as at late times. As seen from Table HH since in this case 
the two non vanishing eigenvalues of the linearization ma- 
trix are of oposite sign, the matter-dominated solution is 
always a saddle equilibrium point of Eqs. (|14c[) . 

Equilibrium points P2, P3, and P4, are associated with 
early time dynamics since, according to Eq. (|16p . v = 1 
is correlated with infinitely large values of the variable <f>. 
Besides, in this case, 



d 6 In V 



-A. 



V = V e 



-\a> 



(28) 



in correspondence with the limit of potential at large 
<j>'s. This confirms that equilibrium points with v = 1, 
indeed correspond to early-times dynamics in the model 
with the cosh-like potential ([5]). The main porperties of 
these equilibrium points (P2 — P4), together with those 
of the matter-dominated phase P\ , have been discussed 
infill]. These can be summarised as follows (the proper- 
ties of the matter-dominated solution Pi have been dis- 
cussed above): 

• The kinetic energy-dominated solution (point Pi in 
Table [T|, is always decelerating. For A < V6 it is 
the past attractor for any phase space trajectory. 

• The scalar field-dominated solution (equilibrium 
point P3), exist whenever A < \/6, and it is infla- 



tionary for A < y/2 (decelerating otherwise). It is 
always a saddle point in the phase space \& defined 
in Eq. |26| . 



• The scaling solution (fixed point P4), exists when- 
ever A > This solution is correlated with 
decelerated expansion of the universe for 7 > 2/3 
(it is inflationary otherwise). The 0- field mimics 
matter with a barotropic index 7^ = 7, while their 
energy densities scale as 



1 



A 2 /3 7 - 1 



(29) 



It is always a saddle point in ^. 

• The late-time dynamics driven by potential |(SJ| is 
correlated with the minimum V m \ n = A at <f> = 
and its neigbourhood -, which, in the phase space 
\& , is depicted by the equilibrium point with v = 
in Table[T] (point P5). This point corresponds to 
the inflationary de Sitter solution 2>H 2 = A. Since 
the real parts of the eigenvalues of the linearization 
matrix for P5 are all negative (see Table [TTJ) , the 
de Sitter solution is always the future attractor for 
any phase space trajectory. 

For values of the parameter A > \fa/2, the de Sitter 
attractor is a spiral fixed point and, whenever a < 24 
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both, the past attractor (point P2), and the future (late- 
time) attractor (point P5) co-exist, meaning that the cor- 
responding orbits in the phase space will be repelled by 
the kinetic energy-dominated solution and will be even- 
tually attracted to the de Sitter solution. 

In this case (A > y/a/2), the linear perturbations of 
the scalar field perform damped oscillations around the 
minimum of the potential at — v = 0, that are 
characterized by cyclic frequency 

For A 2 > a/4, the general solution for the evolution of 
linear perturbations Sxi = (Sx, 8y, Sv) in the neighbour- 
hood of the minimum of the potential (equilibrium point 
rc. m tables Tables U and HJ), can be written as 

8x t = a lie XlT + a l2 e X2T + a l3 e X3T 

= ajl e- 37r + e - 3T/2 (a l2 e MT + a l3 e- MT ) ,(31) 

where Ai, A2, and A3 are the eigenvalues of the lineariza- 
tion matrix around P5. 

The latter oscillatory behaviour with frequency u> is 
what, according to Refs. [Til. H2I [l3| . can be associated 
with cold dark matter because the amplitudes of the per- 
turbations decrease at a rate cx exp (— 3r/2). The corre- 
sponding mean energy density (p^) then dilutes at a rate 
cx a -3 with an effective equation of state (u>(f>) — (see 
Eq. ©). 

The above results are illustrated in Fig.[5]and Fig.[5]for 
a fixed value 7 = 4/3. In particular, the oscillating solu- 
tion is depicted in the latter figure, where it is apparent 
the way the orbits of the ODE (|14c[) . at late times, coil 
around the segment {{x, y, v) = (0, y, 0) : < y < 1} un- 
til, eventually they reach the inflationary de Sitter attrac- 
tor P 5 = (0, 1,0). The oscillating solution arises due to 
the choice of the free parameters in Fig. [5] A = 5, a = 3, 
that obey the constraint A > ^fa/2. 

For self-consistency, we show in Figs. |H] and [7J the 
late-time (point P5 above, associated with the neigh- 
bourhood of the minimum of the potential), and early- 
to-intermediate times local behaviours of the orbits of 
Eqs. (|14c[) , respectively, for small values of the parame- 
ter a (0 < a < 1). 

Due to the choice of the values of the free parame- 
ters (A = 2, a = 0.1), obeying A 2 > a/4 as in the for- 
mer case, since the corresponding equilibrium point is a 
spiral one, the orbits coil around the de Sitter segment 
(x,y,v) = (0,y,0) until, eventually, they reach the de 
Sitter solution 3H = A (Fig. The early times-to- 
intermediate behaviour (Fig. [7]), clearly shows that the 
kinet energy-dominated solution (equilibrium point P2) 
is the past attractor, while the matter-scaling solution 
(point P4) is a saddle equilibrium point. 

We want to emphasize that the above spiral form of 
the orbits of the ODE (|14c[) - see Fig. [6] - is what can 
be properly correlated with CDM behavior, so that, only 
for A > \fa/2 the scalar field component in our model 
may be called as SFDM. 




0.2 0.4 0.6 0.8 
«« 



FIG. 4: Orbits of the autonomous system of ODE (|14c|l with 
function f(v) in Eq. (|22p . and the corresponding flux in time 
r - upper panels -, for the chosen values of the parameters 
(A = 2, a = 10, 7 = 4/3- background radiation), for different 
sets of initial conditions. In the lower panel we show the 
projection of the orbits into the phase plane (x,y). Notice 
that the kinetic energy-dominated solution (a;, y, v) = (1, 0, 1) 
is the past attractor, while the de Sitter solution (0, 1, 0) is the 
future attractor. The radiation-dominated solution (0, 0, v) is 
always a saddle fixed point for any given v. 



V. SUMMARY OF RESULTS 

The relevant results of the study of the fixed points 
of the autonomous system of ODE (| 14c[) . for the cosh- 
like potential V(4>) in Eq. ([5]) - corresponding to f(v) in 
Eq. (|22|) -, can be summarized as follows (see the former 
section). 

• The past (early-time) attractor in the phase space 
4", for A < 6, is the kinetic energy-dominated (de- 
celerated) solution 3P 2 = cj) 2 /2 (see Figf4]). For 
values A > 6 there does not exist any past attrac- 
tor in \t (see FigfS]). 

• The scalar field-dominated phase 3P 2 = cp 2 /2 + 
V{<t>) (existing only for A 2 < 6, and inflation- 
ary whenever A < \/2), and the scaling solution 
fi /fi 7 = l/[(A 2 /37)-l] (exists whenever A 2 > 37, 
decelerating if 7 > 2/3), are always saddle equilib- 
rium points. 

• The (decelerating) matter-dominated solution 
3P 2 = p~f exists for any value v in , 3>. It always 
represents a saddle fixed point (see Fig. [4j . 

• The future attractor is always the inflationary de 
Sitter solution 3P 2 = A. At late times, for A > 
\fa/2, the orbits of the ODE (|14c|l coil around the 
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FIG. 5: Different perspectives of the local late-time behaviour 
of the orbits of the autonomous system of ODE (|14c|) . with 
function f(v) in Eq. (|22[) , for different values of the parameters 
A = 5, and a = 3 (7 = 4/3- background radiation) - upper 
panels. In the lower panel we show the projection of the orbits 
into the phase plane (x,v). At late times the phase space 
trajectories coil around the de Sitter segment {(0,j/,0) : < 
y < 1} until, eventually, they reach the inflationary de Sitter 
attractor (0, 1,0). 



FIG. 6: Local behaviour of the orbits of the autonomous 
system of ODE (114c [I with function f(v) in Eq. (l25jl . in the 
neightbourhood of the point P5 associated with late-time cos- 
mic dynamics, and the corresponding flux in time r - upper 
panels -, for the chosen values of the parameters: A = 2, 
a = 0.1 (7 = 4/3- background radiation), for different sets of 
initial conditions. In the lower panel we show the projection of 
the orbits into the phase planes (x,v) and (x,v), respectively. 



de Sitter segment {(0, y, 0) : < y < 1} until, even- 
tually, these reach the late-time de Sitter attractor 
3H 2 = A (see Fig. p. The later (late-time) be- 
haviour can be associated with oscillating solutions 
that reproduce CDM in the model. 



In the next section we will discuss the relevance of the 
above results to the cosmology of the model under inves- 
tigation. 



VI. DISCUSSION AND CONCLUSIONS 

The fact that the scaling solution (point P4 in Tablc[T|, 
the matter-dominated phase (point Pi), and the de Sit- 
ter inflationary solution (P5) are equilibrium points of 
Eqs. (|14cj) for the potential (O, is not a happy coinci- 
dence but a profound principle of nature, inherent in the 
SFDM model with cosh-like potential under investigation 
in this paper. 

In terms of the dynamical systems analysis this means, 
in particular, that independent on the initial conditions 
chosen, there is a chance for the cosmological model to 
transit by each of the mentioned phases, spending some 




FIG. 7: Early times-to-intermediate behaviour of the orbits 
of the autonomous system of ODE (|14c[) with function f(v) 
in Eq. (125(1 . for the chosen values of the parameters: A = 
2, a = 0.1 (7 = 4/3- background radiation), for different 
sets of initial conditions. In the right-hand panel we show 
the projection of the orbits into the phase plane (x,v). It 
is evident that the kinetic energy-dominated solution (point 
P2 in Tab(T]| is the past attractor, while the scaling solution 
(point P4 = (0.81, 0.71, 1)) is a saddle equilibrium point. 



time in the neighbourhood of each one of them. 1 More 



The neccessary requirement for the occurrence of an equilibrium 
point x' = y' = v' = 0, means that the universe may keep evolv- 
ing in the neighbourhood of this point for a quite long time, 
depending on the nature of the fixed point: if the equilibrium 
point is a saddle in the phase space, it is correlated with a tran- 



9 



precisely, independent on the initial conditions, the cos- 
mological evolution might drive the universe from one 
phase to another to, inevitably, end up in the de Sit- 
ter stage (point P5) being the late-time attractor in the 
phase space of the model (5'). 

The cosmic evolution driven by the 0-field after the 
end of inflation, most probably starts in the neighbour- 
hood of P4 where the scalar field tracks matter (radia- 
tion). Then, at early-to-intermediate times it approaches 
to the matter-dominated (decelerated) stage Pi where 
most structure was formed until, finally, the evolution 
end ups in the inflationary de Sitter stage P5. 

The fact that, for appropriate values of the free pa- 
rameters (A > y / a/2), the latter solution is a spiral equi- 
librium point means that, at late times, while the phase 
space trajectories approach to the de Sitter phase, the 
scalar field perturbations perform damped oscillations 

sient stage of the cosmic evolution, but if it is an attractor, it may 
only represent either the starting point of the evolution — past 
or early-time attractor or its end point - future or late-time 



that mimic CDM. Therefore the SFDM model with a 
cosh-like potential could also be a nice scenario to ad- 
dress the united description of the dark matter and of 
the dark energy in a single field [20j. 
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